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I. INTRODUCTION 



Entanglement, together with its apphcations, is one of the main features of quantum in- 
formation theory Q, [3] ■ It is a resource for new communication and computation algorithms. 

A pure state n = oi a. quantum system establishes quantum correlations between 

its subsystems, entangling them with each other. As a general rule, the more mixed (in the 
sense of majorization) the reduced density matrix v:^ = Ttb tt is, the stronger will be its 
entanglement with the other parts. In bipartite quantum system the entanglement is the 
same for either part, and we may speak of the entanglement between both subsystems. In 
addition, if one part is 2- dimensional, the orbits of the reduced density operators under local 
unitary transformations depend on one parameter only. 

The problem of characterising entanglement becomes more difficult when the total system 
is in a general (i. e., mixed) state. There are now quantum as well as classical correlations. 
Their distinction depends on the task in question and is, hence, not unique. Therefore, 
generally, one has to choose between several entanglement measures Is, 4]. Among them, the 
certainly most important one is the entanglement of formation Eq>{p), discovered by Bennett 
et al. [Sj, expressing the asymptotic number of ebits (maximally entangled qubit pairs) 
needed to prepare a given bipartite state p by local operations and classical communication 
(LOCC). Let $ denote a trace preserving positive map from one quantum system into itself 
or into another one, and denote by S'$(p) the von Neumann entropy of the output $(p), 
given the input state p. Then we have 

£;^(p) = min^p,S($(vr,)) (1) 

where the minimum is taken over all possible convex i^Pj = l,Pj > 0) decompositions of 
the state p into pure states 



P 



^PjTTj, TTj pure, i.e., TTj = (2) 



Let us call this quantity entanglement entropy of $ or ^-entanglement for short. This 
provides the entanglement of formation, if $ is specified in Eq. ([T]) to be one of the partial 
traces, Tr^ or Tr^, of a bipartite quantum system. In other words, the entanglement of 
formation is the ^-entanglement with $ = Tr^ or $ = Tr^. The construction above 
preserves the symmetry between both parts of a bi-partite quantum system observed in the 
pure state case. 



A further example for the appearance of the global optimization problem Eq. ([T]) is the 
HSW theorem of Holevo, Schumacher, and Westmoreland . It gives the one-shot or 

product state classical capacity x{^) of ^ channel $ by first subtracting E^{p) from S^{p) 
and then maximizing this Holevo quantity x*{p) over all input density operators: 

xUp) = smp))~EM (3) 

X>s> = max x*M 
p 

Closed formulas for the entanglement of formation, i.e., analytic solutions to the global 
optimization problem Eq. ([1]) are only known for certain classes of highly symmetric states 
[s], ^, for the ^-entanglement of a 3-dimensional diagonal channel |lO|] and for the exceptional 
case of a pair of qubits. In this case of a 2 x 2 system one knows a complete analytic 
for the entanglement of formation. It has been obtained first for rank two states j^. 



'ormula 
and 



later generalized to all 2-qubit states by Wootters 
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Wootters expressed E^{p) in terms of another entanglement measure C$(p), called con- 



currence m 



Generally, one can replace the von Neumann entropy S in Eq. ([T]) by any other unitary 
invariant, preferably concave, function, say G, on state spaces. Substituting G{^{n)) for 
S'($(7r)) in Eq. ([T]) one obtains another entanglement measure attached to positive and 
trace preserving maps. The concurrence is a measure of this kind: Let $ map the states of 
a quantum system into those of a 1-qubit quantum system, i.e., a map of output rank 2. 



Then the ^-concurrence C$ is defined by using G{p) = 2y/det p. To get the concurrence of 
bipartite a 2 x n-system one sets $ = TrR . The concurrence appeared to be an interesting 
entanglement measure in its own right [3]. Many authors, e.g. jl^, Q, 16], have obtained 
bounds for the concurrence of general bipartite systems. 

We may now state the aim of the present paper as follows: We study C$ and for 
general 1-qubit trace preserving positive maps $. We also exemplify in Section HVl D how to 
transform our results to rank two density operators of a 2 x n quantum system. 

In section [Til we explain important properties of roofs and describe, for a positive and trace 
preserving map $ from any quantum system into a 1-qubit system, the relation between $- 
concurrence and <l>-entanglement, including entanglement of formation. In Section UTTl we 
provide an explicit expression for the concurrence of general positive (stochastic) 1-qubit 



maps. We found this construction in 



171]. Afterwards we learned that a similar result had 



already been obtained by Hildebrand In this paper we elaborate on those results. 

The Section IIIII contain a streamlined version of the constructions and proofs of Hildebrand 
and our unpublished work. 

Our construction of the concurrence works for all stochastic (trace-preserving positive 
linear) 1-qubit maps, not only for completely positive ones. It is, therefore, suggestive but 
not the topic of the present paper, to ask for applications to the entanglement witness 
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problem 

Section [IV] is devoted to a more detailed study of examples. We present explicit formulas 
and intuitive pictures of the convex roof construction for some important classes. We start 
with bi-stochastic 1-qubit maps (subsection A), followed by a short discussion of 1-qubit 
channels of Kraus length two. Subsection C explores the richness of stochastic maps com- 
muting with rotations about an axis. The last subsection D explains, mainly by example, 
the application of our previous results to more general channels (trace preserving and com- 
pletely positive maps) with 1-qubit output. Section |V] is devoted to the $-entropy for axial 
symmetric stochastic maps. We find several qualitative different phases distinguished by the 
geometric pattern of their roofs. In Section |Vl] we shortly discuss the use of our construction 
at concurrence problems for channels with higher rank. 



II. THE CONVEX ROOF CONSTRUCTION 



Let us elaborate on some details of the solution of the global optimization problem 
Eqs. f|T|2|) by the so-called convex roof construction. Let G be a function on the convex 
set fl of density operators of a finite quantum system. A point p G is a roof point of G if 
there is an extremal convex combination Eq. ([2]) such that 

G{p) = J2p^G{n,). (4) 

Then the convex decomposition p = ^Pji^j with pj > and YlPj — ^ ^ill be referred to 
as G- optimal. Thus, if we knew a G-optimal decomposition of p, we could calculate G{p) 
from the values attained at pure states. A roof point p will be called flat if there exists an 
optimal decomposition Eq. (jlj) where all values G{nj) are mutually equal, i.e., G{p) = G^iij) 
for all j. 

The function G will be called a roof if every density operator p of is a roof point for G. 
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Similary one defines a flat roof as a function G for which every point p is a flat roof point. 

Let gin) be a function defined on the set of pure states. Then G is called a roof extension 
of g if G{p) is a roof and G{7i) = gi^n) for all pure vr. On the other hand, if Gconv is 
a convex extension of g from the pure states to all states then Gconv ^ G for every roof 
extension G. The assertion can immediately be seen from Eq. (jl]) and the very definition of 
convex functions (Jensen's inequality). Since the supremum of any set of convex functions 
is convex again, there is a largest convex extension which is, however, not larger than any 
roof extension of a function g. Is this largest convex extension a roof? One knows that the 
answer is "yes" for continuous g. Continuity of g, together with the compactness of the set 
of pure states, guaranties that the largest convex extension of 5^ is a roof and, hence, the 



unique convex roof extension of g 
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2l|: 



Theorem 1. Let giji) he a continuous real-valued function on the set of pure states. There 
exists exactly one function G{p) on which can he characterized uniquely hy each one of 
the following four properties: 

1. G is the unique convex roof extension of g. 

2. G{p) is the solution of the optimization prohlem 

G{p)= inf y^PjOi^^j)- (5) 



3. G{p) is largest convex extension of g 122] . 

4- G is the smallest roof extension of g. 

Furthermore, given p E Q, the function G is convexly linear on the convex hull of all pure 
states TT appearing in optimal decompositions of p. 

Therefore, G provides a foliation of Q into compact leaves such that a) each leaf is the 
convex hull of some pure states and h) G is convexly linear on each leaf. 

If G is not only linear hut even constant on each leaf, it is called a flat roof. 

Item 1 of the theorem justifies to write "min" instead of "inf" in Eqs. (jS]) and ([T]). 
Let us apply the theorem to find out how concurrence and ^-entanglement relate for 
stochastic maps $ from an arbitrary quantum system into a 1-qubit system. Setting (a la 
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Shannon)^ H{xi,X2) = — Xilogaji — X2logX2, one has the following: 

Theorem 2. Let $ a stochastic map into the states of a 1-quhit system. Denoting by 
its ^-entanglement and by C$ its concurrence. The function 



e(x)=//^^,^J, l=x^ + y' (6) 
is strictly convex within — 1 < x < 1. It holds 

EM>acM). (7) 

and this is an equality when p is a flat roof point of C$ . 

To prove this theorem we have to collect three facts: a) For pure states tt we have equality 
in Eq. ([7]) and the value of both sides is the von Neumann entropy of $(7r). Hence, both 
sides are extensions of ^(^(Tr)). b) The right hand side of Eq. ([7]) is convex, see appendix 
A for a proof. The left hand side is a convex roof and, hence, not smaller than any other 
convex extension. This proves the inequality Eq. ([7]). c) If p is a fiat roof point of C$, then 
the same is true for any function of (7$, in particular for ^(C$). Therefore, the left hand 
side, being a convex extension, cannot be larger then the right one and equality holds. 

In the case of the entanglement of formation of a 2-qubit system ($ = Tr^) the concur- 
rence is a flat roof and, hence, equality always holds in Eq. ([7]). This has been proved by 
Wootters by explicitly constructing fiat optimal decompositions for all 2-qubit density 
operators. 

However, already the concurrence of a 2 ® 3 bipartite system or of a general 1-qubit 



channel is not a fiat roof. Eq. (JTj) together with the Fuchs-Graaf inequality ( 23| , see also 



2J]) for 1-qubit states S{p) < 2 (log 2) v^det p provides then the estimate 

acM) < EM < log(2) CM (8) 

for all stochastic maps with 1-qubit output space, i.e., for all stochastic maps of (output) 
rank 2. 



^ Our formulas are valid for arbitrary bases of the logarithm. The basis 2 is used for numerical calculations 
and plots of, e.g., the HSW capacity. 
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III. STOCHASTIC 1-QUBIT MAPS 



The space A^2 of hermitian 2x2 matrices p = ( ) is isomorphic to Minkowski space 
via 

x=(xo,x) p = ^{xol + x-a) (9) 

1 / Xo + X3 Xi+ iX2 

2 \xi - iX2 Xq - X3 

We have det p = j{xl — xl—xl — x|) = |x- x where the dot between 4- vectors denotes the 
Minkowski space inner product and Tr p = xq. Therefore the cone of positive matrices is just 
the forward hght cone and the state space f2 of a qubit, the Bloch ball, is the intersection of 
this cone with the hyperplane V defined by a;o = 1- In this picture mixed states correspond 
to time-like vectors and pure states to light-like vectors, both normalized to Xq = 1. 



A trace-preserving positive linear map $ : A^2 ^ -^2 can be parameterized as 
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= $ Q(a;o/ + x-a)^ = ^ {xqI + {xot + Ax) ■ a) (10) 

where A is a 3x3 matrix and t a 3- vector. 

We consider the quadratic form g on A^2 defined by 

4 

g*(x) = 4(det $(p) - wdet p) = $(x) • $(x) - wx • x = ^ QijXiXj (11) 

i,j=0 

where w is some real parameter. For pure states, i.e., on the boundary of the Bloch ball 
where x ■ x = 0, the form g(x) equals the square of the concurrence C = 2^ydet $(p). 

Furthermore, we denote by Q the hnear map Q : Xi J2 lij^j corresponding to the 
quadratic form q via polarization: 

Qt = Qt-y^m,= { ' ^ ^ J (12) 

where rjij = diag(-|-l, — 1, — 1, — 1). The following two statements are the central result of 
this section: 

Theorem 3. Let the quadratic form q and therefore the matrix Q be positive semi- definite 
and degenerate, i.e., Q > and dimKerQ > 0. IfKeiQ contains a non-zero vector n which 
is space-like or light-like, n ■ n < 0, then q^^^ is a convex roof. Furthermore, this roof is fiat 
if such an n exists with Uq = 0. 



Theorem 4. For every positive trace-preserving map $ there exists a unique value wq for the 
parameter w such that the conditions of Theorem\^ are fulfilled. Therefore, the concurrence 
of an arbitrary stochastic 1-qubit map $ is given by C^{p) = a/^1|^Ip)- 

Let us sketch the proof of Theorems [3] and |H The square root y/q of a positive semi- 
definite form g on a hnear space provides a semi-norm on this space and hence it is convex. 
According to Theorem [1] we need to show that it is also a roof, i.e., there is a fohation of the 
space into leaves such that g^/^ is linear on each leaf. Let n = (no, n) be a non-zero vector 
in KerQ. Then for all vectors m we have 



g(m + n) = (m + n)Q{m + n) = mQm = g(m). 



(13) 



Let us start with the case where n can be chosen to have no = 0. Then n gives a direction 
in V along which q is constant. Therefore, ^yq is a flat convex roof. 





















c > 




M2 



FIG. 1: The embedding of the Bloch ball into M.2 and its foliation by a flat convex roof. 



Let us now consider the case where Ker Q does not contain a vector n with no = 0. Then 
we have dim Ker Q = \ and this line intersects V in one point which we call n. Every other 
point m in can be connected to the point n by a line lying in V . Then q^l"^ is linear along 
the half-line R'*' 9 s 1-^ sm + (1 — s)n since 



q (sm + (1 — s)n) = (sm + (1 — s)n)Q(sm + (1 — s)n) 
= s^g(m) 



(14) 



This concludes the proof of Theorem [3l Our proof of Theorem H] presented in 17|] used the 



elegant argument following jis], Il9 | 



Gorini-Sudarshan parametrization 26|| of stochastic maps. Here we give a shorter and more 




FIG. 2: The foliation of the Bloch ball in the case uq ^ 0. 

We will consider the flow of the signature of the quadratic form q = qo — wrj as function of 
w G R. It is clear that for sufficiently large w we have sgng = sgn(— r^) = (+ + H — ) whereas 

for large enough negative w we have sgng = sgn{ri) = (H ). A signature change can 

only occur at one of the real roots Wi of det Q = det(Qo ~ '^'7) = 0. The "Minkowski 
metric" rj is regular and rj = ri~^. Therefore the Wi are the real eigenvalues of tiQq since 
det Q = (det?7) det{r]Qo — wl). 

Positivity of $ implies 5'*(x) > for all x with x ■ x > 0. This is just the assumption 
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of Yakubovich's S-lemma from the theory of quadratic forms (see |l8l . Il9l . l27l | which ensures 
the existence of a non-negative value w such that is at least positive semidefinite, g? > 0. 
Then it is clear that all four eigenvalues of rjQo are real and that Wi > w > W2 > w^, > w^: 
There must be at least one signature change above or at w and at least 3 signature changes 
below or at w. More signature changes are impossible since we have at most four real 
roots. There is (up to degeneracies) only one possible pattern of signature changes and q is 
positive and degenerate, sgng = (+, ..,0), precisely at w = wi and w = W2- It is positive 
definite for wi > w > W2 if wi W2- In the case wi ^ W2 let ni,n2 be the corresponding 
vectors in KerQ^j^. Then xviQqUi = Winf and n2Qon2 = Furthermore, no nonzero 

vector can be both in KerQ^i and KerQ^j since r] is non-degenerate. So, niQoni > W2n1 
and n2Qon2 > Winl (since Qw^^ ^ > 0), providing {wi — ^2)112 < and {wi — ^2)11^ > 0. 
Therefore, KerQ is time-like at Wi and space-like at W2- 

In the degenerate case wi = W2, Ker Q is at least two-dimensional. In this case, let rii, 112 
be two orthogonal (in the Euclidean sense) vectors from KeiQ. Then rii and 112 can not 
both be time-like (since there is only one time-like direction). 

This proofs the claim of Theorem HI existence of a suitable Wq. It is given by W2, the 
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second largest eigenvalue of r/Q* • 

IV. EXPLICIT EXAMPLES 

Let us demonstrate our construction on some examples. From here on we will sometimes 
denote the coordinates xi,X2,X3 of state space Eq. ([9]) as x,y and z. 

A. Bistochastic maps or unital channels 

Bistochastic maps preserve the center of the Bloch ball. We have t = and the Bloch ball 
is pinched by A = diag(Ai, A2, A3). This includes the depolarising channel p 1— > pp+{l —p)^I 
where A = diag(p, p,p) and also the phase-damping channel where A = diag(p,p, 1). We 
get w = max(Af , A2, A|) and 



The Holevo quantity x$(p) (see Eq. [3]) is a concave function. Since the channel is symmetric 
under all 3 reflections Xi i— > — Xj, it must take its maximum, the HSW capacity 




(15) 



which is flat in one direction since one of the terms in the sum vanishes. 



Nevertheless, this case includes channels of all Kraus lengths between 1 and 4. 
Since the roof is fiat, the entanglement entropy is given by 




(16) 



X$ = maxxlip) 



(17) 



p 



I — I 



at the origin of the Bloch ball, p = ^I. This reproduces the well-known 





1 - Vw 
2 



(18) 



B. Channels of Kraus length 2 



A channel has Kraus length two if it can be represented as 



$(p) = A^pA + B^pB 



(19) 
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The concurrence of such channels has already been studied in |29| using a quite different 
approach. According to [s^, unitary transformations can bring such a channel to the form 



A = diag(cosM, cosu, costi cosu) 
t = (0, 0, sin-usinf ) 



(20) 
(21) 



which corresponds to A = [cos | cos |]/ + [sin | sin |]cr^, B = [cos | sin |]o"a; — ^[sin | cos |]cry 
and we can assume cosm > cosf. Then we find for the concurrence w = cos^m and 



C%{p) = ?/^(cos^('u) — cos^(f)) + (2;cosMsint> — cos f sin -u)^ 



(22) 



which is positive semi-definite and independent of x, so we have again a flat roof. All 
channels which arise from a bipartite 2x2 system with rank-2 input states via restriction of 
the partial trace to the support space of the input state are of length 2 and have therefore 



a fiat roof, in accordance with Wootters' celebrated result 
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C. Axial symmetric channels 

Every positive trace-preserving linear map commuting with rotations about the Xs-axis 
is (modulo unitary transformations) of the form 

^^^^^ /«Xoo + (l-7)xn Pxo, \ ^23^ 

y (3xio 7x11 + (1 - a)xoo J 

with real non-negative parameters a, P, 7. The Bloch ball is pinched by A = diag(/3, P,a + 
7 — 1) and then shifted along the a;3-axis hj t = (0, 0, a — 7). 
This family includes many standard channels. Besides the 

• phase-damping channel (length 2, unital) for a = 7 = 1 and 

• the depolarizing channel (length 4, unital) for a = 7, (3 = 2a — 1 
which we already considered, we also find 

• the amplitude-damping channel (length 2, non-unital) for 7 = 1, jS"^ = a. 

Positivity of $ demands 

< a,7 < 1, (24) 
< PL. = 1 + 2a7 - a - 7 + 2^a{l - a)j{l - 7). (25) 
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The first inequality guarantees that north and south pole of the Bloch ball are not mapped 
to the outside, the the second one describes the hmit when the eUipsoid touches the sphere 
at a circle. 




FIG. 3: A map which is at the boundary of the set of positive maps. 

The stronger condition of complete positivity of $ evaluates to 

< a7 (26) 
For the concurrence we have found the explicit expression 

Cl(X) = 4(det - w det{X)) (27) 

with 

w = max(/?^/?^) (28) 
where = 1 + 2a-f - a - - 2^/^^^^ayy(T^^. (29) 

In the case /3 > /3c we have a flat roof whose leaves are in planes perpendicular to the 2;-axis. 




FIG. 4: Leaves of the concurrence oX (5 > f5c- 
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In the other case we have a one- dimensional KerQ generated by n = {1,0,0, zq) with 
Zo = ^^^"^ 7)+Vq(i q) ^ r^Y^^ jg j^Q^ rjnj^g leaves are straight hnes meeting at the 

point Zq on the 2;-axis outside the Bloch baU: 



FIG. 5: Leaves of the concurrence at P < (3c- 

At the bifurcation point (3 = f3c the concurrence is linear everywhere on the Bloch ball 
(and therefore every decomposition is optimal): 

Cf,=f,^ (p) = (v/«(l-a) - ^7(1- 7)) z + + ^7(1 - 7) (30) 

The special case of the amplitude-damping channel a = 7 = 1 and therefore [3 = (3c = 
Pmax belongs to this degenerate situation with 

Cad{p) = (1 + z)^a{l-a) (31) 

Since this channel has length 2, this result is also a special case of eq. (!22|) for u = —v with 
a = cos^ u. The concave Holevo quantity must take its maximum for states on the z-axis 
where we get 

Xad(^) = V (^") + ^ (1 - - ^ ((1 + (32) 

The equation ^ = can be solved only numerically. The resulting capacity is plotted in 
Fig.l 
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FIG. 6: HSW capacity for the amplitude-damping channel as function of the channel parameter a. 



Similar results can be found in 



3l|. 



Let us finally mention that we have no explanation for the striking similarity between 
eqs. f l2S]) and fl2^ . They differ only by the sign of the square root. So, P^^x derived from 
Fig. [3] and this (3^ of the roof bifurcation are roots of the same quadratic equation. 



D. The 2 X n bipartite quantum system 

Here we consider 2 x n systems Ti. = 7i"^®7i^, dimTi"^ = 2. Let 'H2 be any 2-dimensional 
subspace of H and V a unitary mapping of H"^ onto Ti.2- Then 

P^Hp)=V{Ttbp)V^ (33) 

is a 1-qubit channel for all density operators p supported by 7^2- The eigenvalues of 
and of = Tr^ p are the same. Hence, by Eq. (fTTj) and by Theorem HI we are allowed to 
write 

C{pf = 4(det p'^-w det p) (34) 

for all density operators p with support in H2 and with a unique w = w(T[2)- Notice that 
this representation does not depend on the choice of the unitary V in Eq. (!33|) . However, w 
depends on the 2-dimensional subspace TC2- 
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As an illustrating example we choose n = 4 and consider Ti^ as a 2-qubit system. Then 
7i becomes a 3-qubit system, 7i = ®Ti.^ <^H'^, and the partial trace Tr^ from Eq. fl33l) is 
identified with Tibc- An interesting subspace is generated by the W and GHS state vectors 
given by \W) = S-i/^qooi) + |010) + |100)) and \GHS) = 2-i/2(|ooo) + |111)). Defining 
the unitary V in Eq. dSSD by V|0) = \GHZ) and = $ can be computed to be the 

1- qubit map 

This an axial symmetric channel and we can read off w = 1/6, therefore, 

8 4 

C^(p) = g^oo + ^11 + -^^ooxu (36) 
For p supported in our subspace this is equivalent to 

dpf = \{GHZ\p\GHZY + {W\p\Wf + \{GHZ\p\GHZ) {W\p\W) (37) 

After this quite explicit example we return to the more general case of Eq. flM|) . We 
rewrite the 2x2 determinants in Eq. (IMl) by the help of the characteristic equation in terms 
of traces: 

C{pf = 2[(Trp^)2 - Tr((p^)2)] - 2w[{Tr pf - Tr(p^)] (38) 
Polarization of this quadratic form provides (compare Eq. ffTTl) ) the bilinear form 

g^(pi,P2) = 2(1 - w7)(Trpi)(Trp2) + 2 [«;(Trpip2) - (Trpfp^)] (39) 

defined for all pairs of Hermitian operators on Ti. If pi and p2 are supported by the same 

2- dimensional subspace 7-^2, and if w is correctly chosen, then is positive semi-definite 
and degenerate on that subspace. Hence, if C(pi) = 0, then also qw{pi,P2) = for all p2 
supported by H2- In particular, if pi = tti is a separable pure state and p2 a state, we get 

1 - Tr nfp^ = w{l - Tr TTipa) (40) 

It holds TTi = TT^ (g) 7rf , as tti is assumed separable. 

If there is a second pure separable state, say tt2, supported by H2, one gets 

1 - Trntn^ = w{l - (TrTr^^^) (Trvrf vrf )) (41) 
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Thus, in this particular case, the number w is determined by the transition probabihties 
Tr 7rf ''^TT^''^ = \{^i''^\4'2'^)\'^ between the marginal states of tti and 712. One observes that 
w can vary between and 1 already for subspaces generated by two separable vectors. This is 
a nice illustration of Theorem [31 The operator — tti belongs to Ker Q, and the concurrence 
remains constant along the intersection of the Bloch ball carried by 7^2 with every real line 
of the form p + t{TT2 — tti). 



V. ENTANGLEMENT ENTROPY FOR AXIAL SYMMETRIC STOCHASTIC 1- 
QUBIT MAPS 

In this chapter we study the entanglement entropy defined in Eq. for the axially 
symmetric map Eq. (l23l) in more detail, using Theorem [2] and numerical methods. 

Our aim is an understanding of the structure of the foliation of the Bloch ball provided 
by the convex roof construction. This foliation encodes the optimal decompositions Eq. ([2]) 
for all states. The foliation changes with the channel parameters. In most of the 7) 
parameter space all states have an optimal decomposition into two pure states. In a small 
region of the parameter space we find optimal decompositions of length 3. We characterize 
the bifurcation structure of this "phase transition" and its position in parameter space. 



There exist quite a 



qubit channels, e.g. 
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ot numerical and analytical work about the HSW capacity of 1- 



33, 



34| where the optimal decomposition of the optimal state is 



considered. In contrast, we consider the optimal decomposition of all states. 



A. Some degenerate channels 

a. a = 7 In this case the channel is unital and has therefore a flat convex roof for the 
concurrence. We have /3max = 1 and = (2a — 1)^, so we find w = max((2a — 1)^, /?^). 
The concurrence C, and hence too, are constant either (in case of {2a — 1)^ > on 
concentric cylinders around the 2;-axis -E^ = Eq,{x'^ + y"^) or on planes perpendicular to the 
z-axis = Ei^{z). 

b. a + 7 = 1 In this case the range of the channel is degenerate, being a 2-dimensional 
ellipse orthogonal to the 2;-axis. Furthermore, Pc = and therefore w = (3"^. We get again a 
fiat roof. C<^{z) and hence E'^, too, are constant on planes perpendicular to the z-axis. 
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B. The general case (a — 7) (a + 7 — 1) 7^ 

We did extensive numerical studies of the global minimization problem of the entangle- 
ment entropy Eq. ([T]), guided by and compared to analytic studies of special cases. The 
following overall picture emerged: There are 3 different phases. For fixed values of a and 
7, we have at large values of /3 a phase (phase I) where the entanglement depends only on 
z. By decreasing /5, we reach phase II where a cone with apex at the north pole appears. 
States in the cone have optimal decompositions of length 3. The opening angle of the cone 
decreases and for small enough /3 we reach phase III, where again all optimal decompositions 
have length 2. 



Phase la Phase lb Phase II Phase II Phase II Phase III 




/3max>/3>/3c /3c>/3>/3i /3i>/3>/32 /3i>/3>/32 /3i > /5 > /52 /32 > > 
FIG. 7: Leaves of the foliation of the entanglement entropy. The z axis points upwards. 



Of course we have a flat entanglement roof as long as we have a flat concurrence roof 
(phase la). But the flat phase for the entanglement extends to even lower values of /3, where 
the concurrence is not longer flat (phase lb)! 

For phase III let us remark that the leaves form cones with their apex on the 2;-axis 
outside the Bloch ball. But different to the Phase II of the concurrence (compare Fig. [5]) 
they do not intersect at the same point on the 2;-axis. 

The above picture and the equations for /3i and /32 below are valid in the case 

(a-7)(a + 7-l) >0. (42) 

For the opposite case, turn the pictures upside down [z —z) and exchange a ^ 7 in the 
equations below for /3i and /32. 

The bifurcation points /5i and ^2 between the 3 phases can be calculated analytically. 
Let s(cos(0)) denote the entropy ^(^(Tr)) for the pure state vr = (sin(0), 0, cos(0)). Then 
the bifurcation point /?i can be found by comparing the competing decompositions Ex = 
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|s(l) + fs(cos(0)) with^a = cos(0)). We expand Ei(0) -^2(0) = ^(a, /3, 7)0^+0(0^) 

and get Pi as the root of g{a, P, 7) = 0. 

Using the abbreviations x = 2a — l,y = 2'j — 1 v/e find 



Pi =777 — — ^ — Y\ : — TTT (^^ + xy + (a;^ - l)l/arctanh(x)- 

2[x + (x^ — Ij arctanh(xj ^^g-^ 

— x^) arctanh(x)(x^ — xy^ — (x^ — arctanh(x)) j 

Analogously, we obtain /?2 by comparing the decompositions Ei = s{l) + 

i-co&((^) ^^_^^ _ ^(•(.Qgj^^^^ around (f) = n: 



2 ^. (1 + x) log(l - + (1 - x) log(l + y)-{l + x) log(l + x) - (1 - x) log(l - x) 



2(log(l-y)-log(l + y)) 

(44) 



C. Phase diagram 



The following figure shows the phases in the (3, 7-plane for a = 0.8. The upper boundary is 
given by the positivity condition, Eq. (1251) . The boundary between phases la (entanglement 
and concurrence have fiat roofs) and lb (only entanglement has fiat roof) is given by Eq. ( l29l) . 
Phase II is bounded by Eqs. P3|) and 
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ca 




FIG. 8: Phase diagram in the (7, /3)-plane for a = 0.8. 

The phase II region where length 3 optimal decompositions exist as well as the phase 
lb are quite small but they exist everywhere outside the degenerate points where either 
a + 7 = lora = 7. 

D. One-shot (HSW) capacity 

The Holevo quantity will take its maximum for a state on the z-axis. Its numerical 
calculation is highly simplified by taking the foliation structure into account. We show in 
Figure [9] the P dependence of this maximum, i.e., the HSW capacity, for fixed values of a 
and 7. The maps are positive for (3 < Pmax, completely positive for (3 < Pep- The values Pi 
and P2 indicated in the figure separate the phases I, II and III. In phase III the capacity is 
independent of p. 
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0.9 




I 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P 



FIG. 9: The HSW capacity as function of the channel parameter /3 at a = 0.8, 7 = 0.4. The inset 
shows the small region where the transition between phases I, II and III takes place. 



VI. CONCURRENCE FOR CHANNELS WITH HIGHER INPUT OR OUTPUT 
RANK 

Our method provides a complete solution for the concurrence of trace-preserving positive 
maps of input and output rank 2. How could one possibly overcome the input rank two 
(or 1-qubit map) restriction? The following problem may be of interest: Assume Xl^'i^j 
is an optimal decomposition for the concurrence of a 2 x system. Every pair ttj, tt^ of 
different pure states is supported by a 2-dimensional Hilbert space Tijk- Hence there is 
a number wjk = w{Hjk) defining the concurrence for density operators supported by Tijk 
according to Eq. flM|l . Which restrictions on the set of all wjk arise from the optimality of 
the decomposition? 



Another issue is the generalization to higher output ranks. Rungta et al [35j proposed to 
replace the determinant det p by the second elementary symmetric function of the eigenval- 
ues, C$(7r) = 2A/e2($(7r)). While the square root of 62 is concave, one might find a value 
for w making the expression 

2v/e2($(p))-w;e2(p) (45) 
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a convex extension of 2e2($(vr))^/^, vr pure. In these cases, the expression Eq. (1451) is a 
lower bound for the ^-concurrence. An example is the diagonal map Dm in any dimension 
m which cancels the off-diagonal elements. Denoting the matrix elements of p by Xjk-, this 
recipe results in 

Cn{p) > 2(5^ (46) 
j<k 

Another example is the following family of indecomposable Choi maps of a 3 x 3 system: 



^ Xqq + pX22 —Xoi —Xq2 ^ 



1 + P 



(47) 



— Xio Xii + /iXoo —Xi2 
\^ -X21 X22 + pXn J 

is trace-preserving, positive and indecomposable for p > 1. The map $[1] is extremal 
in the set of positive maps. Here our recipe provides the bound 

4u 

C^ipf > (Y^^ [(^00 + xu + X22Y + (/i - 1) (koil^ + |xo2r + ki2r)] , (48) 

a positive semi-definite quadratic form in the matrix entries. In the special case yU = 1 our 
recipe provides an exact though highly degenerate answer: $[1] maps all pure states of the 
3x3 system to mixed states with the same ^-concurrence and therefore the ^-concurrence 
is constant everywhere, C$(p) = 1. 

VII. CONCLUSIONS 

We have explained a way to get concurrences of stochastic 1-qubit maps and of rank two 
states in 2 X n quantum systems. The methods is attractive by its simplicity, providing a 
arge area of applications. The new methods is different from that of Wootters [l^ and of 



36| which is based on conjugations. 

The advantage of the new methods is its applicability to roofs which are not flat. Only a 
small subset of the stochastic 1-qubit maps actually has a ^-concurrence which is a flat roof. 
For a general 1-qubit map the concurrence is real linear on each member of a unique bundle 
of straight lines crossing the Bloch ball. The bundle consists either of parallel lines or the 
lines meet at a pure state, or they meet at a point outside the Bloch ball. Furthermore, C$ 
turns out to be the restriction of a Hilbert semi-norm to the state space. 
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For the special case of an axial symmetric 1-qubit channel we presented a throughout 
study of the ^-entanglement. Here the structure of the optimal decomposition of states can 
be quite different depending on the channel parameters. There is a phase where all optimal 
decompositions have length 2 and are flat, a phase where states with optimal decompositions 
of length 3 exist, forming a cone in the foliation of the Bloch ball, and a phase where all 
optimal decompositions are of length 2 but not fiat. We found explicit formulas for the 
bifurcation points which separate the phases. Interestingly, there exists a region in the space 
of 1-qubit maps where the ^-entanglement is fiat despite the fact that the ^-concurrence is 
not flat. 

Our method of finding optimal decompositions for the concurrence works perfectly for 
rank two density operators only. For higher rank states it provides lower bounds. It is a 
challenge to find an algorithm, if existing, which combines the merits of this approach and 
the conjugation based one. 



APPENDIX A 



The function defined in Eq. ([6]) 



ax)=H[^^-±y], l = x' + y' (Al) 



2 ' 2 

is defined on — 1 < x < 1 and does not depend on the sign of x. It is strictly convex since 



^ + ^ + ] > (A3) 



y^ \ 3 ' 5 ' 7 

Therefore, ^ is the supremum of a family of functions ax + b. Inserting a convex function 
C{p) with values — 1 < C < 1 represents ^(C) by a supremum of convex functions aC + b. 
This proves the convexity of ^{C{p)) as a function of p. 
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